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Preface 


This  thesis  was  originally  conceived  as  a study  of 
scattering  by  circular  cylinders  with  end  caps,  apertures, 
and  internal  cavities  with  a wire  passing  through  the 
cavities.  After  some  research  and  discussion  with  scholars 
knowledgeable  in  this  area  of  study,  it  was  decided  to 
accept  the  suggestions  of  Dr.  Chalmers  M.  Butler  and  limit 
the  analysis  to  scattering  by  two  concentric  finite  cir- 
cular cylinders. 

I am  deeply  grateful  to  my  thesis  advisor  Dr.  William 
A.  Davis  who  made  the  degree  of  success  that  was  attained 
possible  by  spending  many  hours  explaining  the  concepts 
involved  and  how  to  apply  them.  Finally,  I owe  special 
thanks  to  my  family  who  have  sacrificed  tremendously 
that  I might  complete  this  work  and  in  particular  to  my 
wife  Mary  Jo  for  typing  and  editing  this  manuscript. 
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Abstract 


This  thesis  deals  with  electromagnetic  scattering 
from  two  concentric  finite  length  cylinders.  The  scat- 
tering bodies  are  assumed  to  be  perfect  electric  con- 
ductors and  the  incident  field  is  assumed  to  be  a plane 
wave  whose  magnetic  field  is  transverse  to  the  longi- 
tudinal axis  of  the  scatterers.  Electromagnetic  field 
equations  are  derived  in  terms  of  vector  potentials. 
Integral  equations  are  developed  for  the  surface  currents 
induced  on  the  cylinders.  A numerical  solution  is  form- 
ulated in  terms  of  the  method  of  moments  for  the  unknown 
surface  currents  and  the  components  of  the  scattered 
electric  field.  Results  are  presented  for  scatterers 
of  different  dimensions.  Comparisons  are  made  between 
scattering  from  a finite  hollow  cylinder  and  scattering 
from  two  concentric  finite  cylinders. 
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I.  Introduction 


Survivability  in  a nuclear  war  environment  is  a facet 
of  weapon  system  design  that  must  be  considered  if  the 
weapon  system  is  to  be  effective  in  such  an  environment. 

In  addition  to  other  effects,  one  must  consider  the  effect 
of  the  electromagnetic  pulse  (EMP)  generated  by  a nuclear 
detonation.  EMP  may  induce  large  currents  in  a weapon 
system  which  damage  or  destroy  critical  electronic  compon- 
ents, possibly  rendering  the  system  ineffective.  Thus, 
one  needs  to  be  able  to  predict  these  currents  and  design 
the  system  to  minimize  their  effects.  There  are  three 
approaches  that  could  be  used  to  obtain  EMP  data.  One 
could  detonate  a nuclear  device  and  measure  the  induced 
currents.  This  is  obviously  impractical  for  environmental 
and  political  considerations.  Secondly,  one  could  build  a 
simulator  which  produces  electronically  a wave  similar  to 
an  EMP.  The  Air  Force  Weapons  Laboratory  has  an  EMP  sim- 
ulator which  is  used  to  obtain  test  data.  Unfortunately, 
it  is  impractical  to  test  all  systems  in  this  manner.  The 
final  approach  to  predicting  EMP  effects  is  to  formulate 
the  problem  as  an  electromagnetic  scattering  problem. 

The  electromagnetic  scattering  approach  is  pleasing 
in  that  one  could  determine  the  EMP  effects  on  different 
designs  without  actually  building  the  different  systems 
and,  theoretically,  any  system  could  be  evaluated.  However, 
the  number  of  structures  for  which  the  scattering  problem 
can  be  solved  analytically  is  severely  limited.  In  general, 
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analytic  solutions  are  limited  to  small  classes  of  canonical 
structures.  The  advent  of  the  high  speed  digital  computer 
has  permitted  the  analysis  of  electromagnetic  scattering 
by  more  general  structures  through  numerical  approximation 
methods.  Although  much  work  has  been  devoted  to  this  approach, 
many  more  structures  need  to  be  simulated  before  an  accurate 
prediction  can  be  made  of  the  EMP  currents  induced  on  com- 
plex structures  such  as  aircraft  or  missies. 

This  thesis  represents  an  effort  to  analyze  the  electro- 
magnetic scattering  from  two  concentric  finite-length  cylinders. 
Davis  (Ref  4)  investigated  electromagnetic  scattering  by  a 
finite  hollow  cylinder.  The  problem  considered  here  will 
be  restricted  in  a manner  similar  to  Davis's  by  considering 
only  perfect  electric  conductors  and  transverse  magnetic 
plane  wave  incidence.  An  effort  is  made  to  present  the  dif- 
ferences between  scattering  from  a finite  hollow  cylinder 
and  the  scattering  from  the  same  cylinder  with  an  inner 
cylinder,  or  wire,  whose  radius  is  much  smaller  than  the 
radius  of  the  outer  cylinder. 

The  information  in  this  thesis  is  grouped  into  chapters 
based  on  the  traditional  subdivisions  and  the  theoretical 
and  mathematical  concepts  that  must  be  considered  in  solving 
this  type  of  problem.  Chapter  II  contains  the  theoretical 
development.  Electric  and  magnetic  field  equations  are 
developed  in  terms  of  vector  potentials  instead  of  the  elec- 
tromagnetic field  quantities  themselves.  This  permits  one 
to  assume  a lower  degree  of  smoothness  on  electromagnetic 
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fields  while  maintaining  a rigorous  development.  The  field 
equations  are  then  manipulated  to  give  integral  equations 


for  the  equivalent  currents  induced  on  a scattering  surface 
by  an  incident  field.  The  electric  field  integral  equation 
is  finally  specialized  to  scattering  by  both  the  concentric 
cylinders  and  the  hollow  cylinder. 

Chapter  III  describes  how  one  might  solve  the  equations 
developed  in  Chapter  II  on  a digital  computer.  The  method 
of  moments  is  described  briefly  and  then  applied  to  the 
scattering  current  equations.  This  results  in  a system  of 
simultaneous  linear  algebraic  equations  whose  solution  rep- 
resents an  approximation  of  the  induced  currents.  Some  man- 
ipulations are  shown  which  facilitate  solving  this  system 
of  equations  on  a digital  computer.  Finally,  digitized  ex- 
pressions for  the  total  electric  field  are  developed. 

Chapter  IV  presents  plots  which  illustrate  the  validity 
of  the  results  and  the  differences  between  the  scattering 
from  a hollow  cylinder  and  the  scattering  from  the  concentric 
cylinders.  Comparisons  are  made  between  Harrington's  solu- 
tion for  scattering  from  an  infinite  cylinder  (Ref  7)  and 
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scattering  from  finite  cylindrical  structures.  Plots  are 
included  for  cylinders  whose  dimensions  simulate  structures 

r 

of  interest  to  the  Air  Force. 

Chapter  V presents  conclusions  that  can  be  drawn  about 
scattering  from  concentric  cylinders  relative  to  scattering 
from  hollow  cylinders.  Finally,  recommendations  are  made 
for  a series  of  problems  which  would  extend  this  work  and 
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provide  data  that  would  be  useful  in  the  design  of  weapon 
systems . 

Three  references  were  particularly  useful  in  the  devel- 
opment of  thi  thesis.  The  dissertation  by  Davis  (Ref  4)  and 
Harrington's  book  (Ref  7)  on  electromagnetic  field  theory  were 
used  extensively  in  the  theoretical  development.  Harrington's 
book  (Ref  6)  on  the  method  of  moments  was  used  for  analysis 
of  how  to  digitize  the  theoretical  equations  for  computer 
implementation.  Finally,  the  series  of  reports  published 
by  the  Air  Force  Weapons  Laboratory  as  Electromagnetic  Pulse 
Interaction  Notes  were  used  to  gain  insight  into  the  magni- 
tude of  the  EMP  problem  and  the  geometries  that  have  been 
investigated  previously. 
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II . Theoretical  Development 


Integral  expressions  for  the  electric  E and  magnetic 
H fields  that  exist  in  an  electromagnetic  scattering  prob- 
lem can  be  derived  in  terms  of  a magnetic  vector  potential 
A and  an  electric  vector  potential  F (Ref  7*77)  • In  gen- 
eral, fields  can  be  expressed  as 

C = c£  + £*  a) 

where  C represents  the  total  field,  C1  represents  the  in- 
cident field  or  the  field  that  would  be  present  if  there 
were  no  scatterer,  and  Cs  represents  the  field  due  to  the 
scatterer.  The  scattered  fields  in  electromagnetic  scattering 
problems  are  due  to  equivalent  electric  and  magnetic  cur- 
rents induced  on  the  scattering  body  by  the  incident 
field.  When  appropriate  boundary  conditions  are  imposed, 
determining  the  total  electromagnetic  fields  reduces  to 
solving  an  integral  equation  for  the  unknown  equivalent 
currents  on  the  scattering  body.  In  this  chapter,  vector 
potentials  and  Maxwell's  equations  are  used  to  develop 
integral  equations  for  the  equivalent  currents  on  a finite 
hollow  cylinder  and  the  equivalent  currents  on  a finite 
cylinder  and  wire,  the  bodies  being  perfect  electric  con- 
ductors and  the  wire  being  on  the  longitudinal  axis  of 
the  cylinder. 

Vector  Potential  Equations 

The  bases  of  all  electromagnetic  field  problems  are 
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Maxwell's  equations.  For  simple  matter,  that  is,  matter 
that  is  homogenous,  linear,  isotropic,  and  time  invariant, 
Maxwell's  equations  can  be  written  in  differential  form 
(ej65t  time  variation  suppressed)  as 


V*  H 

s jojtE  +■  f 

(2) 

7*  E 

- -j  t oju.  JT  - M 

(3) 

= 4/e 

(4) 

V * H 

- rm/M 

(5) 

where  E and  H are  electric  and  magnetic  field  intensities, 

J and  M are  electric  and  magnetic  current  densities,  and^ 
and /m  are  the  electric  and  magnetic  charge  densities. 
Taking  the  divergence  of  Eqs  (2)  and  (3)*  substituting  into 
Eqs  (4)  and  (5),  and  noting  that  the  divergence  of  a curl 
is  zero,  yields  the  continuity  equations 


V • j = 

(6) 

y • s - 

(7) 

An  expression  for  the  electromagnetic  fields  in  terms 
of  the  magnetic  vector  potential  A can  be  derived  by  assum- 
ing there  are  no  magnetic  current  sources,  M = 0,  and  taking 
the  divergence  of  Eq  (3).  This  yields  V-H'O  . Pierce 
(Ref  10s221)  states  that  for  any  vector  whose  divergence 
is  identically  zero  everywhere,  there  exists  a vector 
whose  curl  equals  the  original  vector.  Therefore,  let 


6 


H = 7 * fl  (8) 

where  A is  defined  as  a magnetic  vector  potential.  Sub- 
stituting for  H in  Eq  (3)  yields 

7*  (t  *■  jcJM  R ) - O (9) 

For  any  vector  whose  curl  is  identically  zero  everywhere, 
there  exists  a scalar  whose  gradient  equals  the  vector 
(Ref  10 i 221).  Therefore,  let 

E - }cOM  A = “ V (10) 

where  V is  defined  as  an  electric  scalar  potential.  Sub- 
stituting for  E in  Eq  (2)  yields 

7 * 7*  R - j*°£(-  R - 7/)  + J (11) 

Helmholtz  (Ref  10»221)  showed  that  a vector  is  completely 
specified  to  within  a constant  when  both  its  curl  and  di- 
vergence are  given.  To  complete  the  specifications  of  A, 
let  the  divergence  of  A be 

7R  • -jw«_V  (12) 

Substituting  Eq  (12)  into  Eq  (11)  and  using  the  vector 
identity 

V*  7«  « = 7(7-  fl)  - 7*fl  (13) 

results  in 

♦ h'fl  * - J (14) 
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where  k =td>*£'.  Substituting  Eqs  (8)  and  (14)  into  Eq  (2) 
yields 

•jco l£  - It  R +V(v  fi)  (15) 


If  there  are  only  magnetic  current  sources,  J = 0,  one 
can  write  by  the  principle  of  duality  (Ref  7:99) 


vaf  + h'F  = -ft 

(16) 

and 

E --  -V*F 

(17) 

+V(V-F) 

(18) 

where  F is  defined  as  an  electric  vector  potential.  If 
both  electric  and  magnetic  current  sources  are  present,  one 
can  write  by  the  principle  of  superposition 

£--&[?  R)J  - V»  F (19) 

'v*  (20) 

Thus,  once  A and  F are  known,  the  electromagnetic  fields 
E and  H can  be  determined. 

Scattering  Equations 

Equations  (19)  and  (20)  represent  the  total  electric 
and  magnetic  fields  in  simple  matter  obtained  from  independent 
vector  potentials.  If  a scattering  body  is  placed  in  the 
medium,  the  solution  for  the  fields  becomes  somewhat  more 
complicated  since  the  vector  potentials  are  composed  of 
incident  and  scattered  components.  One  must  calculate  the 
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scattered  vector  potentials  prior  to  solving  for  the  elec- 
tromagnetic fields. 

The  vector  potentials  A and  F due  to  finite  indepen- 
dent sources  in  simple  matter,  enclosed  in  a volume  V,  at 
any  point  r in  all  space  are  (Ref  7:100) 


flCFJ  Jv' 

y 

(21) 

Hr)  rjfi  $(?,?')  Jv>' 

(22) 

y 


where  the  prime  notation  indicates  "with  respect  to  the  source 
coordinates".  Note  that  use  has  been  made  of  the  three  dimen- 
sional free  space  Green’s  function 


-j  M r - r'l 

<L 

^nlr-r'l 


(23) 


The  scattered  components  of  the  vector  potentials  for  a 
scattering  body  can  also  be  shown  to  be  given  by  the  forms 
Eqs  (21)  and  (22).  This  may  be  shown  by  considering  the 
scattering  body  in  simple  matter  as  depicted  in  Figure  1. 
All  independent  sources  J are  in  volume  V which  is  bounded 
by  the  surface  of  the  scatterer  3 and  the  surface  which 
approaches  infinity.  The  unit  normal  vector  n at  S points 
into  V. 
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Figure  1.  Scattering  Geometry 


To  determine  the  scattered  field  components,  consider  the 
identity  (Ref  11:300) 


* f*  ■ v$  ♦ #}  * l y*  R(0]j 

in  ^ 


where  the  observation  point  r is  not  on  the  closed  regular 
surface  ^V,  the  boundary  on  V formed  by  3 and  3^.  It  is 
assumed  that  the  field  A meets  the  Sommerfeld  radiation 
conditions  (Ref  2:20),  is  continuous  with  a continuous 
first  derivative  on  S,  and  has  a continuous  second  derivative 
in  V in  the  sense  of  distributions  (Ref  1:114).  Hence,  the 
surface  integral  over  S.  is  zero  and  Eq  (24)  may  be  written 


+ RL*)V'$  + [v-  R(F)]J  iw  r 

+ [*•*(*)]▼$  it  (25) 

Using  Eq  (14)  and  the  expression  (Ref  8:340) 

v'$  <26) 

where  f(r  - r*)  is  the  Dirac  Delta  Function,  a general- 
ized function  or  distribution  as  described  by  Wylie 
(Ref  12»317).  one  can  write 

fl(r)  z j J J Jv  «■  || f *'  r[v'»  ff(r')]  * [*'  •/!(?')]?'§ 

V S 

fl(F')]  - V'f  xRCrOjJ  it'  (2?) 

where  the  reference  system  has  been  changed  by  interchanging 
the  prime  and  unprime  notation. 

One  needs  and  to  determine  E(r) 

and  H(r ) . The  identities  (see  Appendix) 

7*^7'$  * [S' »#<?')]  J,*  r w'*  ItHr')  J*’ 

3 S 

- vfv-jf  A'  X fl  (28) 

*"  s 

vlv-  ffi'}  x [7'*  P-^‘)J  J*'|  =0  (20) 

can  be  used  to  show  that 

+ ${[&*  *(7  *Rt»)) 'V’Jl'f 

» j 

- J - 7')*$  J (30) 
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and 


and  V x ff(r)  ; V*J$  f /»■'  + ? * j § A' X 
V s 

‘•JUffei  *[r'Lv'W‘)J  j.' 

--L-  v{v- Jmf 

J ^ 5 

Thus,  the  expression  for  the  total  electric  field  is 

=■  A*/ $ ^ Jw'  * * (v  j f j 

* v v' 

+ f {$ «'  * £*'» /?(*-;]  - /j'|  [?'• 

-v'}  * p'  * h?‘)]]j3‘  f $ p * [v'*R(f‘j]-vj v'{  J*' 

& 

- £ [*'-R(p)]  (*'•?',)  *'$  J-' 

5 

which  may  be  rewritten  as 

jwt  Etr)  = +rl'V'fi,'Lt)]  * k*  Jf  [v'xflc?')]  J *• 

s 

♦ 7 w'  y v'/V'-  jj(P)]4*‘  + f ft' xfr'xR  (?')]-?' I v‘f  J*' 

s s'-  J 

* 7*1  1**1  S'  x ii  (i‘)  J<' 

s 

where  A1(r)  is  the  magnetic  vector  potential  due  to  the 
independent  source  J and  A(r')  is  the  total  magnetic 
vector  potential  at  S.  The  preceding  equation  may  be 
written  in  simpler  form  by  using  the  unit  dyadic  I (Ref 
9*569) 
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J s 

r7'r']-ff(r)Js-  is‘  04) 

S 

Using  the  same  notation,  one  may  write  the  magnetic  field 
equation  as 

H(f)  r V*  ftLt)  + &<?)] 

s 

* [a'J  + vv]'  f jt  *[kti*v'r']-ffv)  **'  (35) 

s 

If  the  only  independent  source  is  a magnetic  current 
source,  then  by  duality  E(r)  and  H(r)  are 

f (T)  = -Vf  f*Cr)  J J<' 


and 

jUf/i  =[kkI  +w]'F\p)  vf** 

I 

+?'?]  &’)  +[k'I+w]'(§  &'*[?' *hf)Us' 

s 

(37) 

• 

When  both  independent  magnetic  and  electric  sources  are 
present,  the  electromagnetic  field  equations  are 

iw*ik  R‘‘f) 

♦ jyj  - V7]-  h *■  * ’■*  »<'-)  * t Iff  "**l  HJ-' 

♦ /fe£  jVf . vV]-fi(?) -v’*r  (?)}<!•' 

s ^ 

(38) 

and 

} I 
| ( 

f 1 

A»Cf)  * «-rj^  [|ff  + V7j  • TV) 

■ l 

♦ *V?J  { «*-  .{pi  * ??]■  HI?)  ♦ 

s 

? 

■ 

I 

{.  . 

♦ $ V’  Kjy'x  fhr')  [k%I  ♦y'y'J  -Rr')  | Jj' 

(39) 

tvc 

i ’* 

l 

* i 

t 

> 
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Rewriting  Eqs  (38)  and  (39).  one  obtains 


£(»•  E'lf) 


+ Vx’l  $ »'  * E (.?’)  As' 

s 


(40) 


and  H(f)  = H\?)  - jVf  + VV]  xEi?') 

s 

* *»$$*■  M) 


where  E1(r)  and  Hx(r)  are  the  incident  electromagnetic 
fields  or  fields  due  to  the  independent  sources. 

Since  the  equivalent  surface  currents  on  a scatterer 


are  given  by  (Ref  71IO6) 

j*  = *■  x a ('•> 

(42) 

and  fY  = - x £ 

(43) 

the  electromagnetic  fields  in  terms  of  these  equivalent 
currents  are 

Eif.)  = EV)  «•  [>!  ■•'VVj.Jf  Js(?o 

5 

<!•'  (44) 

s 


and  H(f)  * S‘Cf)  Jt' 

s 

+ V*f$  JUr-)  J«' 

s 


(45) 
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If  the  scatterer  is  assumed  to  be  a perfect  electric  con- 
ductor, these  equations  simplify  to 


and 


E (r>  * [a1  i +vv]f  jit  <J 

5 


s 


s' 


(46) 

(47) 


since  Ms  is  zero  on  a perfect  electric  conductor. 

It  should  be  noted  that  the  development  in  this  chapter 
used  vector  potentials  to  reduce  the  restrictions  on  the 
differentiability  of  the  E and  H fields;  the  differen- 
tiability of  the  integrals  in  Eqs  (46)  and  (47)  must  still 
be  considered  when  the  observation  point  is  permitted  to 
be  on  the  scattering  surface.  This  point  will  be  addressed 
in  the  last  section  of  this  chapter. 

Cylinder  Plus  Wire  Equations 

Figure  2 depicts  the  geometry  of  the  scattering  prob- 
lem considered  in  this  thesis.  The  cylinder  and  wire 
(infinitesimally  thin,  finite,  cylindrical  surfaces  at 
/ = a and  / = b respectively)  are  assumed  to  be  perfect 
electric  conductors.  The  wire  is  assumed  to  have  a neg- 
ligible radius  relative  to  the  radius  of  the  cylinder  and 
the  wavelength  of  the  incident  field.  The  incident  field 
is  a plane  wave  assumed  to  be  a transverse  magnetic  field 
with  its  propagation  vector  in  the  x-z  plane  at  an  angled^ 
with  respect  to  the  z axis  as  shown  in  the  inset  of  Figure 
2.  This  ensures  that  there  will  be  a component  of  the 
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incident  E field  along  the  wire. 

In  general,  the  incident  fields  and  currents  for  the 
structures  of  interest  can  be  expanded  in  Fourier  series. 
For  plane  waves,  the  incident  E and  H fields  can  be  rep- 
resented as  (Ref  4tl 7) 


* os 


(48) 


and 


H- 


= £ /> h;  Kh  ^ 


i"4>  kt  e-Oi  4% 


(49) 


» i - oO 


where  the  coefficients  are  a function  of  only  and  the 
variation  with  respect  to  z and  <f>  (coordinate  variation 
in  standard  cylindrical  coordinates)  are  shown  explicitly. 
The  surface  currents  on  the  structures  may  similarly  be 
expanded  as 


X (i‘) 


-L 

!,»•  «0 


* 


(*) 


(50) 


It  will  be  assumed  that  there  is  no  significant  <f> 
variation  in  the  surface  current  on  the  wire  since  its 
radius  is  negligible  relative  to  wavelength.  Therefore, 
the  wire  surface  current  may  be  represented  by 
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(51) 
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Since  the  z component  of  the  incident  E field  is  the  com- 
ponent which  induces  a current  on  the  wire,  one  may  assume 
that  the  incident  E field  has  only  a z component  which  may 
be  expanded  as  (Ref  4 i 17) 

w — \fi6  -ft*  cos 

- ZL  nn  &i)  £ £ (52) 

» : — tO 

where  EQ  is  the  amplitude  of  the  incident  E field  at  the 
origin  of  the  coordinate  system  and  the  Jn  are  Bessel  func- 
tions of  the  first  kind.  The  current  on  the  cylinder  con- 
tains both  terms  of  the  expansion  given  in  Eq  (50) • 

Davis  (Ref  4),  in  his  analysis  of  scattering  by  a 
finite  hollow  cylinder,  showed  that  currents  of  harmonic  n 
are  associated  with  only  incident  fields  of  the  same  har- 
monic and  that  by  symmetry  only  non-negative  values  of  n 
need  be  considered.  Davis's  results  and  the  fact  that 
the  objective  of  this  thesis  is  to  determine  the  effect 
of  the  wire  on  scattering  by  finite  cylinders  permits  one 
to  use  only  the  zero  harmonic  term  of  the  expansions  for 
the  incident  field  and  cylinder  currents  since  the  zero 
harmonic  behavior  will  be  the  dominant  behavior,  Therefore, 
the  equations  considered  in  this  thesis  for  the  currents 
on  the  wire  and  cylinder  and  the  incident  field  are  re- 
spectively 

rw it')  = £ J.w d')  (5D 
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(53) 


?vj--£  t;  <*•)  * i j;  <*•; 

« “o 

and  _ ; A e cos 

r * s/a/  «►.  S/A<  *i.  ) ^ (54) 

Substituting  Eqs  (51) » (53) » and  (54)  into  Eq  (46)  and 
realizing  that  vector  potentials  only  have  components  along 
the  direction  of  their  sources,  one  obtains 

_ _ __  . ;4r*  cos 

ti(f)  ~ lc  **  +L  C. 

+ (kx  * Jl)  j § (*')  As'  (55) 

£ 


where  d is  the  derivative  with  respect  to  z and  J rep- 
z z 

resents  the  total  zero  harmonic  z-directed  surface  current 
on  the  structures.  The  effect  of  the  ff  component  of  the 
cylinder  current  is  zero  since  the  partial  derivative  of 
(z')  with  respect  to  / is  zero. 

Noting  that  the  tangential  E fields  must  be  zero  at 
the  boundaries  of  perfect  electric  conductors,  one  can 
write  the  equations  for  the  currents  on  the  cylinder  and 
wire  as 

-j ;kz.  cot 

"j t*>£  T0(  ko~  sitf  ^i,)  c.  a,*  — 
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(56) 


and 


■j  f„ 


HU  ». 


) £ 


■jk* 


COS 


JV;  J*' 

+ { </a' j 


(57) 


where  a and  b refer  to  the  radii  of  the  cylinder  and  wire 

respectively  and  Su,  and  S_  refer  to  the  surfaces  of  the  wire 

w c 

and  cylinder  respectively.  The  equation  for  the  correspond- 
ing current  on  a hollow  cylinder  is 

T0  {ka.  c,u  ) e SiN«o- 


(J«fc  ♦ ' 


(58) 


Kernel  Evaluation 

The  development  of  the  integral  equations  for  the 
equivalent  surface  currents  on  a finite  cylinder  and  wire 
was  restricted  by  not  permitting  the  observation  point  to 
be  on  the  surfaces  of  the  scatterers.  The  complete  solution 
for  the  fields  requires  that  the  observation  point  be  per- 
mitted to  be  on  the  surfaces.  This  presents  an  added 
difficulty  since  the  kernel  of  the  integrals  in  Eqs  (56), 
(57),  and  (58)  possess  singularities  for  some  observation 
points.  If  the  order  of  integration  and  differentiation  is 
interchanged  in  these  equations,  as  will  be  done  in  Chapter 
III,  the  kernels  have  nonintegrable  singularities  in  the 
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Cauchy  sense  for  some  observation  points.  In  this  section 
a method  is  discussed  for  extracting  these  singularities, 
thus  partitioning  the  kernel  into  a term  which  can  be  treated 
analytically  and  a residue  which  can  be  integrated  numerically. 

The  method  of  development  in  this  thesis  consisted  of 
using  vector  potentials  to  derive  integral  equations  which 
contained  terms  of  the  form 

T(a/ j &-v)  r Ck%  *0  J t-fj  Ter)  d*'  (59) 

s 

Terms  of  the  same  form  would  have  evolved  if  the  develop- 
ment had  been  based  on  E and  H (Ref  4).  The  vector  potential 
method  was  used  since  its  validity  did  not  require  the  elec- 
tromagnetic fields  to  be  twice  differentiable.  However,  this 
method  does  not  reduce  the  degree  of  smoothness  required  on 
the  kernel  of  Eq  (59).  Both  methods  require  the  same  degree 
of  smoothness  of  the  kernel  and  require  the  extraction  of 
the  nonintegrable  singularities  if  T (/,/’  t-*')  is  to  be 
determined. 

Davis  (Ref  4*8)  treated  the  nonintegrable  singularity 
problem  through  a limiting  process  by  letting  the  observa- 
tion point  r approach  a surface  point  r' . The  surface  was 
divided  into  two  regions,  a small  patch  that  contained  r', 
and  the  rest  of  the  surface.  That  portion  which  excluded 
the  small  patch  contained  no  singularities  and  required  no 
special  handling  as  r approached  r' . However,  as  r approached 
r*  in  the  integral  over  the  small  patch,  it  was  necessary 
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to  take  the  Hadamard  principal  value  (finite  part)  of  the 
integral  (Ref  4i89).  The  use  of  the  finite  part  of  this 
integral  is  justified  since  the  part  approaching  infinity 
cancels  with  the  part  of  the  integral  over  the  curve  bound- 
ing the  small  patch  which  also  approaches  infinity.  Hence 
Eq  (59)  may  be  written  as 

£ U*  + J*  ) f f-*0  JlO  (60) 

s 


where  the  double  bar  through  the  integral  sign  is  Davis's 
notation  for  Hadamard  principal  value  and  the  order  of  in- 
tegration and  differentiation  has  been  interchanged.  The 
assumption  that  this  interchange  is  valid  is  not  justified 
here;  however,  it  is  u~eu  to  illustrate  a form  that  occurs 
when  a valid  interchange  is  made  in  Chapter  III.  Equation 
(60)  contains  a residue  and  a singularity  integrable  in  the 
Hadamard  sense . 

One  may  rewrite  Eq  (59)  as 


•o  JT 

tu,/;,*)  *(r+j;)  Jruo  j 


~TT 
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A<  «/j 


(61) 


where 


ft  = I F - r - ) - 7cw)v  + s'  * f - ■* 


(62) 


and  the  bulk  current  I(Z')  is 


x w) - 2*/ C -1**0 


Ji'U  W* 
, i «■  i > 


(63) 


and  <x  is  the  angle  between  the  vector  to  the  observation 
point  and  the  vector  to  the  source.  The  bar  through  the 
integral  sign  implies  Cauchy  principal  value.  Consider 
the  integral 


>r 

K = l Sl - ^ 


(64) 


This  integral  gives  a singular  kernel  about  R = 0.  A 
technique  for  extracting  the  singularity  (Ref  11)  is  to 
make  a transformation  of  the  integration  in  cat  to  an  in- 
tegration in  y,  where  y = sin  «w/2  , and  then  add 

and  subtract  the  Jacobian  of  the  transformation  from  the 
kernel  of  the  integral.  Operating  on  Eq  (64)  in  the  manner 
just  described  yields 


K-J 


C.  - Coi  -i 


- V 


J **'ip 


Ip*  -]/^»  + te-s'j*’  ♦ t/-s' 


(65) 


The  first  integral  in  Eq  (65)  is  nonsingular  and  can 
be  integrated  numerically,  and  the  second  integral  can  be 
integrated  in  closed  form  to  obtain  (Ref  5:150) 
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^i2***/^  [ * Vce-D*  +(r  */r  ) 
- I*  ( - ztf?'  + ]/(e'*'J'  +-  \ f ♦/')*'  ) jf 


For  the  geometry  considered  in  this  thesis,  one  must 
evaluate  where  iTF'  approaches  zero.  To  facilitate 
this  evaluation  on  a digital  computer,  one  rewrites  as 


* n 


— ■ ) 

^ L \ -■•)*♦  (/vj-/ 


- I*  ( * - 


^ ■ -,7 


Expanding  the  natural  logarithms  in  Taylor  series  yields 


. _J f 2= 

k,  --  n>  £ 


Thus,  one  obtains 


f -i'fcR 


♦ — , A* 


/*i  ■*  l/(f  J 
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or 
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8”  * [(.-.■/  +(,./■)]¥  zi" 


(70) 


when  fir  approaches  zero. 

Applying  the  results  of  the  kernel  evaluation,  the 
equations  for  scattering  currents  on  the  cylinder  and  wire  are 


_j  V*  cos  &i 

~~  j Ca3  L E0  ( k CL  S IV  ) £•  SIM  ^ 2 
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where  a and  b are  the  radii  of  the  cylinder  and  wire  re- 
c w 

spectively  and  I (z')  and  I (z1)  are  the  total  currents 
on  the  cylinder  and  wire  respectively.  The  expression 
for  the  zero  harmonic  current  on  a hollow  cylinder  is 
obtained  by  setting  Iw(z')  equal  to  zero,  yielding 


~j  E0  Jc  (ha-  sji*  •*.  ) 


£ Sj A/  &■  s. 


e. j 4 

•81*  K 


+ H>  (*,<)]  W 


(73) 


In  the  case  of  scattering  by  the  wire  and  cylinder 
there  are  two  equations  and  two  unknowns,  Ic(z')  and  Iw(z'). 
The  case  of  scattering  by  the  hollow  cylinder  has  one 
equation  and  one  unknown  Ic(z').  Since  these  equations 
cannot  be  solved  analytically,  an  alternate  method  must 
be  used  to  obtain  an  approximate  answer.  Chapter  III 
discusses  one  method  of  obtaining  approximate  solutions 
for  equations  of  this  type. 
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III.  Numerical  Development 


The  solution  of  the  scattering  problem  considered  by 
this  thesis  requires  determining  unknowns  that  appear  in 
the  kernels  of  singular  integral  equations.  One  method  for 
solving  linear  equations  of  this  type  is  the  method  of 
moments.  This  process  requires  expanding  the  unknown 
function  in  terms  of  a set  of  expansion  or  basis  functions. 
Each  term  in  the  set  is  operated  on  in  the  same  manner 
as  was  the  unknown  function.  A matrix  equation  involving 
the  coefficients  of  the  expansion  can  be  found  by  taking 
the  moment  of  the  forcing  function  and  of  the  operation 
on  each  expansion  term  with  respect  to  the  members  of  a 
set  of  weighting  functions.  This  chapter  includes  a brief 
description  of  method  of  moments  and  develops  digitized 
forms  of  Eqs  (71)  through  (73)  which  can  be  solved  on  a 
digital  computer.  The  material  presented  in  this  chapter 
is  based  on  the  work  of  Harrington  (Ref  6)  and  Davis  (Ref  4). 
method  of  Moments 

Consider  the  linear  operator  equation 

L u = f (74) 


where  L is  a linear  operator,  u is  an  unknown  function, 
and  f is  a known  forcing  function.  If  the  solution  to  Eq  (74) 
is  unique  and  exists  for  all  f,  one  can  write 

,:L"  f 
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(75) 


where  L~^  is  the  inverse  operator  and  it  if,  assumed  that 
the  domain  of  L and  u and  the  range  of  f are  defined  over 
appropriate  spaces. 

To  solve  Eq  (75)  by  using  the  method  of  moments, 
expand  u in  terms  of  a complete  set  of  expansion  functions 

u - u 


(76) 


where  the  an  are  the  expansion  coefficients.  In  computer 
applications,  the  summation  is  truncated  and  Eq  (76)  becomes 

A/ 

U.  = £ **  «■  (77) 

»z  i 

where  e is  the  truncation  error.  Equation  (7^)  can  be 
written  as 

A/  „ . 

21  L u*  = r - Le,  (78) 


If  one  defines  the  interproduct  of  two  functions  h^  and 
h2  as 

^1,  = J A ,6»)  Ax 


(79) 


and  defines  a suitable  set  of  weighting  functions 
w , m = • • • »M,  Eq  (79)  can  be  expressed  as 


(80) 
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If  the  scalars  ^wm«  Le^>  are  set  equal  to  zero  for  each 
m,  Eq  (80)  reduces  to 

ZL  o-»  <^j  Luu  ) - 


or  in  matrix  notation 


where 


I m\  m ~ ^ wi»  j Lw.) 

L = f > 


If  M = N and  1 is  nonsingular,  the  solution  of  this 
matrix  equation  is 


and  the  unknown  u is  approximated  by 


u = [u„]  [avj 


where 


[M  : ,%] 
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Therefore,  the  approximate  solution  may  be  written  as 


W'  [Q 


(88) 


An  appropriate  set  of  expansion  functions  and  weighting 
functions  must  be  chosen  before  Eq  (88)  can  be  applied  to 
Eqs  (71)  through  (73).  Davis  (Ref  3*39)  suggests  that  the 
generalized  spline 

(2‘)  r sin^a-  1 z - Zj)  ?;A  (i‘-  2„)  (89) 


where 


O 


i Jt'-sJ  <-* 
<2,  It'  -t*  I Z * 


(90) 


is  an  advantageous  choice  for  the  expansion  function  in 
problems  of  this  type  and  that 

( *'  “ *-)  •S  Mn  A (9D 

are  appropriate  weighting  functions. 

Digitized  Equations 

Using  the  expansion  function  of  Eq(89),  the  currents 
on  the  wire  and  cylinder  can  be  approximated  by 

N 

jw  (i'3  = 2L  i*w  s/a>  *(*  - $4  (*''0  (92) 

* * * i 
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and 


k (4- 1 t'-zj)  rZA  (*'-?„) 


(93) 


V 

I ‘(e')-L^  *•* 
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L 

where  A = N+l  . 

Substituting  Eqs  (92)  and  (93)  into  the  scattering 
current  equations  derived  in  the  previous  sectionf  one  ob 
tains  equations  of  the  fc~m 

_ a'  A * *1 

- j U3  £ J o ( A *.  S l*i  d;-  ^ c.  S//V  - 
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where  the  Hadaraard  form  has  been  used  as  a result  of  the 

interchange  of  the  integration  in  z'  and  the  operation 

(k^  + d^)  and 
z 

c ; *-*’)*  J £ z±±±-±l>  U +k,  («>) 
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(94) 


(95) 
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Equation  (94)  can  be  integrated  by  parts  twice  with  respect 
to  z’  to  obtain  an  equation  of  the  form 


_ j 4 * * o-i. 

E-  O Ta  C i *-  61*  dj  ) C.  SIN  = 
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where  the  operation 
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w:  I 


(97) 


may  be  taken  analytically  to  give 


) 4 
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Cot  (44) 


(98) 


To  complete  the  method  of  moments,  one  takes  the  inner 
product  of  the  pulse  weighting  functions  and  Eq  (96)  and 
the  similar  equation  that  represents  the  evaluation  of  the 
boundary  condition  at  the  wire.  Operating  on  the  forcing 
functions  yields 


m\A  * 4/t 

f ( a. ) - ~ J j t ( 4 4_  4 / 
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f = - \ jc^e.  Eo  To  (A  b </v  £ 
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Similar  operations  on  the  right  hand  side  of  Eq  (96) 
results  in  a set  of  2M  linear  equations  for  the  2N  unknovm 
current  coefficients.  These  equations  have  the  form 
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where 
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for  the  cylinder  plus  wire  problem.  These  equations  can 
be  written  in  matrix  notation  as 
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The  resultant  solution  for  the  current  coefficients  is 
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if  the  matrix  is  nonsingular  and  invertible.  The  solution 
for  the  hollow  cylinder  problem  is 
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where  1 (a, a)  and  f (a)  are  defined  in  the  same  manner  as 

mn  m 

the  definitions  in  the  cylinder  plus  wire  problem. 

Expressions  for  the  total  electric  field  can  be  written 
in  terms  of  these  unknown  currents.  Consider  the  z compon- 
ent of  the  electric  field. 

= e;  + ev>  e‘  (o  doe) 


where  E (a)  and  E (b)  are  components  of  the  field  scattered 
z z 

by  the  structures  at  radii  a and  b respectively.  The  field 
scattered  by  the  structure  at  radius  a is  given  by 

J (.0  - ( k | $ j*  c/s'  (109) 
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Noting  that  this  equation  is  of  the  same  form  as  was 
the  equation  that  had  to  be  solved  in  determining  the 
currents,  one  can  write 
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function  of  the  currents  or  the  components  due  to  the  cur- 
rent on  each  of  the  N divisions  of  the  structures  in  the  z 
dimension  and  t\  is  the  impedance  of  the  medium.  The  total 
E^(r)  is  given  by  a summation  over  n for  both  structures  or 
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where  the 
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(112) 


and  n = 1,2,3»*«*»N. 

An  expression  for  the  component  of  the  electric 
field  can  be  derived  using  Eq  (4),  the  equation  for  the 
divergence  of  E,  yielding 
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When  the  observation  point  is  on  the  wire,  Eq  (113)  can  be 
written  as 


s \s  "IT  UiM 


since  E is  zero  at  the  surface  and  E*  is  always  zero  for 

£d  * 

the  geometry  considered  in  this  thesis,  lising  Eq  (6), 
the  current  continuity  equation,  one  can  write 
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about  the  radius  j - b.  Integrating  from  _/>  = b-  to  s = b+ 
yields 
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and  Ey,(b)  becomes 
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Assuming  that  E^,  is  zero  inside  the  wire,  b, 
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where  central  differences  about  z = nA  has  been  used  to 
obtain  . Likewise,  E,  for  / = a is 
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or 
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For  observation  points  not  on  either  surface 
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which  may  be  written  in  difference  form  as 
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or 
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The  term 


p Jfi  . o 

^ i*  24 


(124) 


must  be  interpreted  in  some  manner  with  respect  to  j . 
Since  a difference  form  was  used  on  , one 
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interpretation  is  to  use  an  average  with  respect  to  to 
give 


aI‘ 


(*) 


- eKy.iJ 

4A 


(125) 


In  this  chapter,  equations  that  can  be  solved  digitally 
have  been  developed  for  the  unknown  scattering  currents. 
Additional  equations  were  developed  for  the  components  of 
the  total  electric  field  in  terms  of  these  scattering  cur- 
rents and  the  incident  field.  Chapter  IV  describes  some 
of  the  specific  structural  dimensions  that  have  been  sim- 
ulated on  a digital  computer  and  the  results  of  these 
simulations . 
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IV.  Results 


The  objective  of  this  thesis  was  to  determine  the 
effect  of  placing  a second  scatterer  inside  of  a finite 
circular  cylinder.  Two  computer  programs  were  written  to 
provide  data  for  this  analysis.  One  program  computed  the 
zero  harmonic  z-directed  currents  and  the  electric  field 
components  for  a hollow  cylinder.  The  second  program  com- 
puted the  same  quantities  for  a cylinder  of  the  same  di- 
mensions as  the  hollow  cylinder;  however,  a wire  was  placed 
on  its  longitudinal  axis.  The  programs  were  run  for  structures 
of  various  dimensions  in  an  effort  to  provide  data  which 
demonstrates  the  validity  of  the  development  and  programming 
as  well  as  data  which  would  be  useful  in  analyzing  practical 
scattering  problems.  The  output  included  plots  of  the  in- 
cident E field,  the  z-directed  currents  on  the  structures, 
and  the  total  electric  field  components  for  five  different 
z planes. 

Validation  Data 

Although  all  results  were  analyzed  for  effects  that 
differed  from  what  one  would  expect,  two  specific  geometries 
were  simulated  to  check  the  validity  of  this  work.  The  first 
case  simulated  a cylinder  whose  radius,  a,  made  the  product 
ka  equal  to  the  first  zero  of  the  zero  order  Bessel  function 
of  the  first  kind.  The  second  case  simulated  a long  thin 
cylinder. 
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One  can  readily  analyze  the  effect  of  the  cylinder 
by  selecting  ka  to  be  equal  to  the  first  zero  of  Jg(k_f  ) 
and  assuming  normal  incidence.  As  would  be  expected,  there 
was  no  scattering  from  the  hollow  cylinder  since  the  forcing 
function  was  zero  at  the  surface.  The  cylinder  current 
was  zero  and  the  total  field  was  equal  to  the  incident 
field.  Adding  the  wire  to  the  geometry  modified  the  results. 
In  addition  to  the  current  induced  on  the  wire,  a small 
current  was  induced  on  the  cylinder  by  the  field  scattered 
by  the  wire.  The  total  scattered  fields,  composed  of  a 
tangential  component  Ez  and  radial  component  E j>  , were 
non-zero  with  the  wire  present. 

Harrington  (7s233)  analyzed  scattering  by  an  infinite 
cylinder,  assuming  normal  TM  incidence,  and  derived  the 
equation 


T*  - zEo 


(126) 


for  the  zero  harmonic  current  density  induced  on  the  cylin- 
der  where  Hq  (ka)  is  the  zero  order  Hankel  function  of  the 
second  kind.  Figure  3 illustrates  the  magnitude  of  the  cur- 
rent computed  by  the  program  for  a hollow  structure  with  a 
radius  of  0.01  wavelength  and  a length  of  five  wavelengths 
and  the  magnitude  of  the  currents  computed  by  using  Eq(133)- 
The  amplitude  variations  of  the  program  computed  current 
are  due  to  the  end  effects  associated  with  a finite  structure. 
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Figure  3.  Magnitude  of  I for  Infinite 

z 

' - 

' •$ 

and  Long  Finite  Cylinders 
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If  one  takes  the  average  of  the  amplitude  variations, 
there  is  good  agreement  with  Harrington's  current  on  an 
infinite  cylinder.  Figure  4 shows  Harrington's  solution 
for  E of  an  infinite  cylinder  and  the  program  solution  for 

Z 

E in  a plane  mid-way  between  the  ends  of  the  cylinder, 
z 

As  one  can  see,  there  is  no  noticeable  difference  in  the 
fields.  Comparing  these  results  to  the  results  for  shorter 
cylinders,  the  convergence  of  the  program  solutions  toward 
the  infinite  cylinder  solution  was  observed  as  the  length 
increased . 

Practical  Data 

The  sponsors  of  this  thesis  are  interested  in  pre- 
dicting the  scattering  effects  from  complex  structures 
such  as  an  aircraft  with  its  myriad  of  cables,  compart- 
ments, equipment,  and  apertures.  In  an  effort  to  provide 
data  that  might  give  insight  into  scattering  problems  of 
this  type,  three  sets  of  dimensions  were  selected  that 
simulate  physical  structures  being  irradiated  by  a TK 
wave  of  unit  magnitude.  The  wavelength  of  the  incident 
field  was  selected  to  be  in  the  middle  of  the  VHF  band  at 
3 0 meters. 

The  first  structure  simulated  was  a coaxial  trans- 
mission line  whose  outer  radius  was  three-eighths  (3/8) 
inch  and  inner  radius  was  one-eighth  (l/8)  inch.  The 
bulk  current  on  the  outer  conductor  was  three  orders  of 
magnitude  greater  than  the  bulk  current  on  the  inner 
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conductor.  Comparing  the  coaxial  transmission  line  to  a 
hollow  cylinder  (circular  waveguide)  with  the  same  outer 
-radius,  one  finds  that  the  z-component  of  the  electric 
field  inslde--the  coaxial  line  is  three  orders  of  magnitude 
smaller  than  the  corresponding  component  in  the  hollow 
cylinder;  whereas  the  radial  comporierrfe -of  the  electric 
field  inside  the  coaxial  line  is  three  orders  of 'magnitude 
greater  than  the  corresponding  component  in  the  hollow 
cylinder.  Although  the  magnitude  of  the  total  current  on 
the  coaxial  line  was  approximately  equal  to  the  magnitude 
of  the  current  of  the  hollow  cylinder,  the  interior  field 
structures  were  significantly  different. 

The  second  structure  simulated  was  a tunnel,  radius 

eight  feet,  with  a conduit,  radius  two  and  one-half  inches, 

on  its  central  axis.  The  length  of  the  tunnel  was  selected 

to  be  75  meters.  Figure  5 depicts  the  magnitude  of  E for 

z 

the  hollow  tunnel  and  tunnel  plus  conduit  as  a function  of 
radius  in  a plane  approximately  9 meters  into  the  tunnel. 
One  can  see  that  there  is  no  noticeable  difference  in  the 
two  cases  and  that  Ez  is  approximately  zero  inside  the 
tunnels  (radii  less  than  2.4  meters).  Figure  6 depicts 
the  magnitude  of  the  radius  times  the  radial  component  of 
the  electric  field  as  a function  of  radius  for  the  same  z. 
One  can  see  that  the  hollow  tunnel  has  no  appreciable  field 
inside  the  structure;  however,  the  tunnel  with  the  conduit 
along  its  central  axis  has  a significant  radial  electric 
field  in  its  interior. 
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The  third  physical  structure  simulated  was  a large 
aircraft  or  missile  with  a radius  of  approximately  six 
meters  and  a length  of  seventy-five  meters.  Figure  7 
depicts  E for  a z-plane  approximately  19  meters  from  the 
end  of  the  structure.  There  is  no  noticeable  difference 
between  the  hollow  structure  and  the  structure  which  had 
a 0.3  centimeter  radius  wire  on  its  axis.  Figure  8 depicts 
the  magnitude  of  radius  times  the  radial  component  of  the 
electric  field  for  the  hollow  structure  and  the  structure 
with  the  wire  in  the  same  z-plane.  Again,  there  is  a sig- 
nificant radial  component  inside  the  structure  with  the  wire; 
whereas,  the  hollow  structure  radial  component  is  approx- 
mately  zero. 

Ten  additional  simulations  were  made  in  which  the 
parameters  wavelength,  physical  size,  and  the  angle  of  in- 
cidence were  varied.  In  general,  the  results  were  similar 
to  the  results  described  in  previous  paragraphs.  An  ex- 
ception occurred  when  the  length  of  the  concentric  cylinder 
structures  was  selected  to  be  less  than  l/2  wavelength  long 
and  the  radii  of  the  outer  cylinder  and  inner  cylinder  were 
selected  to  be  0.1  and  0.0001  wavelengths  respectively.  Fig- 
ure 9 depicts  Ez  in  a z-plane  approximately  three-sixteenths 
(3/16)  of  a wavelength  from  the  middle  of  the  structures. 

One  structure  is  0.5  wavelength  long  and  the  other  is  0.4 
wavelength  long.  It  can  be  seen  that  there  is  a marked 
difference  in  the  two  curves.  The  magnitudes  of  the  currents 
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on  the  outer  conductors  for  the  two  cases  were  approximately 
equal;  however,  the  magnitude  of  the  current  on  the  inner 
conductor  of  the  shorter  structure  was  approximately  five 
times  the  magnitude  of  the  current  on  the  inner  conductor 
of  the  longer  structure.  These  results  suggest  some  type 
of  TM  modal  structure  for  the  fields  in  the  shorter  structure 
and  a TEM  modal  structure  for  the  fields  in  the  longer 


structure . 


V.  Conclusions  and  Recommendations 


This  thesis  addressed  electromagnetic  scattering  by- 
two  concentric,  finite  length,  circular,  perfect  electric 
conductors.  Total  electromagnetic  field  equations  were 
derived  in  terms  of  vector  potentials.  This  permits  one 
to  lower  the  degree  of  differentiability  required  on  the 
electromagnetic  fields  as  opposed  to  what  is  required  when 
the  total  fields  are  derived  in  terms  of  the  electromagnetic 
field  components.  The  total  field  equations  were  manipulated 
to  yield  integral  equations  for  the  equivalent  surface 
currents  induced  on  the  cylinders  by  the  incident  fields. 

The  method  of  moments  was  used  to  formulate  digitized 
equations  that  could  be  solved  on  a digital  computer  for 
the  unknown  currents.  Digitized  equations  were  derived 
for  the  tangential  and  radial  components  of  the  scattered 
electric  field.  Several  comparisons  were  made  between  the 
scattering  from  a finite  hollow  cylinder  and  the  two 
cylinder  case. 

A conclusion  that  must  be  drawn  from  this  analysis  is 
that  introducing  the  second  conductor  into  the  problem 
significantly  changed  the  scattering  problem.  While  the 
magnitudes  of  the  zero  harmonic  currents  on  the  hollow 
conductor  and  the  outer  conductor  of  the  concentric  case 
were  approximately  equal,  the  effect  of  a small  radius 
inner  conductor  cannot  be  ignored.  The  modal  structure 


of  the  fields  between  the  concentric  conductors  is  sig- 
nificantly different  from  the  modal  structure  of  a hollow 
conductor  of  similar  dimensions.  If  one  desires  to  analyze 
the  fields  between  two  concentric  conductors,  the  formula- 
tion must  include  the  scattering  from  the  inner  cylinder 
and  the  cross-coupling  between  the  two  structures.  Further- 
more, if  one  desires  to  design  a cylindrical  structure  that 
will  not  permit  significant  propagation  of  electromagnetic 
waves  interior  to  the  structure,  consideration  of  waveguide 
modes  alone  is  not  sufficient.  Significant  TEM  modes  can 
exist  when  an  inner  cylinder  is  part  of  the  scattering 
problem . 

The  Air  Force  Weapons  Laboratory's  objective  to  predict 
scattering  from  structures  much  more  complicated  than  the 
geometry  considered  in  this  thesis  leads  one  to  suggest 
a series  of  research  projects  related  to  this  work.  A 
suggestion  for  further  study  is  to  simulate  end  caps  on 
the  cylinders  and  various  shaped  apertures  in  the  side  of 
the  outer  cylinder.  One  could  further  extend  the  analysis 
by  placing  a cavity  or  enclosure  behind  the  aperture  with 
the  inner  cylinder  passing  through  the  cavity.  This  analysis 
could  be  further  extended  to  include  multiple  apertures, 
cavities,  and  inner  conductors.  Finally,  non-symmetrical 
scatterers  should  be  analyzed.  This  line  of  research  would 
lead  to  a much  more  accurate  prediction  of  the  scattering 
effects  associated  with  complex  structures. 
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Appendix 


Two  identies  were  used  in  Chapter  II  which  enables 
one  to  write  the  curl  and  the  gradient  of  the  divergence 
of  the  magnetic  vector  potential  A(r)  in  a convenient  form. 
This  appendix  demonstrates  the  validity  of  these  identities. 
The  restriction  placed  on  the  identities  was  that  the  ob- 
servation point  was  not  on  the  surface  of  the  scatterer, 
r / r* . 

Consider  Eq  (28) 

x (S'  xR)  is'  =•  - ^ S'  x k R 

- f (28) 

Interchanging  the  order  of  integration  and  differentiation 
the  left  side  of  Eq  (28)  in  accordance  with  Leibnitz's 
rule  (Ref  12«313)  yields 

| V*[v’$  * (&*  *R)]  Js‘  (127) 
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Note  that 


V(»'-fi)*o  <129) 

and 

v-(3*-p)r0  (130) 


since  (n*  x A)  and  (n' . A)  are  functions  of  the  prime  co- 
ordinates only.  Therefore,  Eq  (28)  simplifies  to 

k'fl  J»*  + f [(£«$)•  rjv'J  - 

3 5 

N*x  ff  «-  v(v- 1 $ w'xff  Ji'j  (131) 

s 3 

where  Eq  (26)  implies  V'V  $ s k $ for  r / r'  and  V$»-V  § , 
Consider  the  expression 

- V (v  • $ { * R J*'j  (132) 

s 

which  may  be  written  as 

-V  + $ V •(»  * R )]  Js ' (133) 
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# 

after  interchanging  the  integration  and  divergence  operators 

or 

-v  $■  (*•*»)•  vj  j*' 

s 

(134) 

i 

since  V(n*  x A)  is  zero  by  Eq  (13 Q).  Interchanging  the 

order  of  integration  and  differentiation  in  Eq  (133)  per- 

mits  one  to  write  the  expression  as 

-f {(*i  v)b3'*A)  / } Ji' 

(1 35) 

Using  Eq  (129)  one  can  write  Eq  (13*0  as 

$[($'*«) -v] 

(136) 

which  completes  the  proof. 

1 

The  second  identity  used  was 

l I 

I 1 

*■ 

v[v- j"'#*  v‘ (v'-R)  ° 

s 

(29) 

f ! 

i i 

which  may  be  rewritten  as  (Ref  7*450) 

i 

! 

Hi  , 

k\. 

vjV  * J V'*  $ v'(v'-  P ) j v'J  = <3 

(137) 

* .< 
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Expanding  the  left  side  of  Eq  (133)  yields 


* ?'(*'•&)  v'*v'(v'-ff)]  Jr‘  ) ZO  (138) 


or 


? (*'•*)  Jv' j ro  (139) 

V 

since  the  curl  of  a gradient  is  equal  to  zero. 

Interchanging  the  divergence  and  integration  operators 
on  the  left  side  of  Eq  (140)  and  expanding  yields 

V f • (v*  v'§  ) - T v X R )]}  Jv'  --0  (140) 

V 1 

which  is  indeed  equal  to  zero  since  the  curl  of  a gradient 
is  equal  to  zero. 
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surface  currents  and  the  components  of  the  scattered  electric 
field.  Results  are  presented  for  scatterers  of  different 
dimensions.  Comparisons  are  made  between  scattering  from  a 
finite  hollow  cylinder  and  scattering  from  two  concentric 
finite  cylinders. 
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